Efficient asymptotic frame selection for binary black hole spacetimes using asymptotic 

radiation 
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Previous studies iiave demonstrated that gravitational radiation rehably encodes information 
about the natural emission direction of the source (e.g., the orbital plane). In this paper, we demon- 
strate that these orientations can be efSciently estimated by the principal axes of (Lf^jL;,)), an 
average of the action of rotation group generators on the Weyl tensor at asymptotic infinity. Eval- 
uating this average at each time provides the instantaneous emission direction. Further averaging 
across the entire signal yields an average orientation, closely connected to the angular components 
of the Fisher matrix. The latter direction is well-suited to data analysis and parameter estimation 
when the instantaneous emission direction evolves significantly. Finally, in the time domain, the 
average (^L(^a.Lb)) provides fast, invariant diagnostics of waveform quality. 



I. INTRODUCTION 



Ground-based gravitational wave detectors like LIGO 
and Virgo are likely to see many few-stellar-mass black 
hole (BHibinaries formed through isolated [l|, IH and dy- 
namical processes. Additionally, advanced detec- 
tors could see the merger signature of two intermediate- 
mass black holes (each M € [100, 1O"^]M0), perhaps 
formed in dense globular clusters Q. Black hole spin 
will significantly impact all phases of the signal to which 
these ground-based detectors are sensitive: the late-time 
inspiral, merger signal, and (through the final BH spin) 
ringdown. Misaligned spins break symmetry in the or- 
bital plane, generally leading to strong modulations in 
the inspiral [7[, merger, and ringdown. At early times, 
these modulations can be easily estimated through sepa- 
ration of timescales: the emitted inspiral waveforms arc 
often well-described by quasistationary radiation from a 
circular orbit, slowly rotating with time modulated with 
time as the orbital plane precesses 0-[ll ■ In the extreme 
but astrophysically important case of BH-NS binaries, 
the emission direction can process about the total an- 
gular momentum through arbitrarily large angles, lead- 
ing to strong amplitude and phase modulations 0, fl0| . 
At late times, numerical relativity simulations of merg- 
ers can efficiently predict gravitational radiation from a 
variety of spins and orbits. A well-chosen static or time- 
evolving frame makes it easier to compare their outputs 
to one another and to these analytic models. 

Previous studies have demonstrated the gravitational 
radiation carried to infinity encodes information about 
the instantaneous emission directions. For example, 
Schmidt et al. Q track the precession of the orbital plane 
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by selecting a frame that maximizes the amplitude of the 
I = \m\ = 2 modes at each time. In this paper, we pro- 
pose an invariant algebraic method to calculate the nat- 
ural principal axes of asymptotic radiation. This method 
can be applied independently at each time, at each fre- 
quency, or at each binary mass. As with previous studies 
[8| the diagonalizcd time-domain representation may be 
useful in interpreting and modeling precessing signals as 
quasistationary circular orbits. Data analysts, however, 
will benefit most from choosing an optimal frame at each 
binary mass. For example, some signals like BH-NS bi- 
naries have an optimal time-domain emission orientation 
that evolves significantly. In these cases, as well as when 
higher harmonics become significant at high mass, the 
frame relevant to data analysis must suitably average 
over the relevant portion, frequency, and modal content 
of the signal. In this paper, we provide the first algorithm 
to select a suitable orientation for the average signal. 

In Section [n] we describe our algebraic method for de- 
termining an optimal orientation versus time, frequency, 
or mass. In Section [m] we demonstrate our method pro- 
vides a new way to efficiently and accurately determine 
the natural orientation of a quasistationary quasicircular 
inspiral versus time, using both artificially generated qua- 
sistationary sources (i.e., applying a rotation to a non- 
precessing source) and waveforms from real precessing bi- 
naries. We also discuss numerical limitations in our time- 
domain method during the merger and ringdown, where 
multiple harmonics become increasingly important. In 
Section IIVI we demonstrate that a preferred frame ver- 
sus mass arises naturally in data analysis applications, 
as the average of the Fisher matrix. Using a separation 
of timescales, we demonstrate this frame roughly cor- 
responds to a time- and bandpass-weighted average of 
the optimal time-dependent frame. We briefly comment 
on the utility of mass-weighted frames for data analy- 
sis applications, such as mode-decomposed searches for 
spinning binaries [ll| . 
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II. NATURAL WAVEFORM FRAMES 

Gravitational radiation carried away to infinity is well- 
known to encode certain preferred orientations, such as 
the radiated linear and angular momenta. More gen- 
erally, given a one-parameter family of inner products 
{a\b)^ = J a*{t')K^{t,t')b{t')dtdt' generated by K^{t,t') 
and a tensor operator acting on the asymptotic Weyl 
scalar, the orientation average of (V'4|Qoi...a„ "04} defines 
an asymptotic tensor 



{Qai...a„)f — 



_ JdQj dtdt'K^{t,t')^l;i*Qa,...a„^p4 



Jdnj dtdt'K^{t,t')ip4*ip4 



■(1) 



The principal axes of this tensor in turn define preferred 
directions. Depending on the inner product and refer- 
ence tensor, very different orientations can be selected; 
this formula, for example, includes as a special case both 
recoil kicks and radiated angular momentum. 

In this paper wc consider (^L(^aLb}}, the average of 
products of rotation group generators L^-^ For exam- 
ple, adopting a kernel Kt{Ti,T2) = 5{ti — t)5{T2 — t) 
(i.e., using a quantum-mechanics-motivated state nota- 
tion, K = \t) {t\), ail average orientation at each time t 
can be calcualted from the principal axes of 



{L^aLb))t 



J dnip4* {t)L(^aLb)ipA[t) 



Eimm' V'4L''04;„ {im' \L(^aLb) \ Im) 
/dr2|V'4|2 



(2a) 



where in the second line we expand i/'4 = 
E;m V'4/m(t)>'i^^^(6','/>) and perform the angular 
integral. In this case, the average tensor <^L((j-^b))( at 
each time t can be calculated algebraically, from the 
mode amplitudes ■04;^ and the well-known action of SU2 
generators Lfc on rotational eigenstates. The principal 
axes follow by diagonalizing this 3x3 matrix. The same 
construct can be applied to any kernel K. In particular, 
we define the average of LaL^ at any frequency / or over 
all frequencies by 



{L{aLb))f ~ 



JdniP4if)L^aLb~)Mf) 

Jdn\MfW 

/dO/d/|^~4|V[(27r/)45,(/)] 



(2b) 



(2c) 



where Sh is the strain noise power spectral density of a 
gravitational wave detector. The subscript M in the sec- 
ond expression denotes the binary mass, used to connect 



^ Explicit coordinate forms for the generators acting on func- 
tions with spin weight s are provided by [12| , denoted as J(. . As 
the algebra of generators acting on the representation is 
unchanged, only the usual action of Lj, on rotation eigenstates 
is needed to perform the calculations described in this paper. 



the scale-free output of a numerical relativity simulation 
to physical time t (i.e., rMtp4{t) depends only on t/M). 
For context, in gravitational wave data analysis, the inner 
product {a\b) = 2 J a*{f)h{f)/ Sh{f) characterizes the 
natural optimal-filtering-induccd metric for a noise power 
spectrum Sh] the denominator in the final expression 
thus corresponds to the inner product {h+\h^) -\- (h^ \hx) 
where the strain J J i/j^dtdt' oc h+ — ihx is reconstructed 
from tpi by integration. 

The components, eigenspaces, and eigenvalues of this 
tensor characterize the emitted radiation, averaged over 
orientation and the specified quantity ^ (= t, /, M, . . .). 
As components of this matrix are associated with aver- 
ages of mode orders weighted by |V'4p, such as 



Ti-{LaLb) = {1(1 + 1)) 



(3) 
(4) 



the components of (^L(^g^Lb)) identify both the dominant 
emission direction and the dominant emission mode. In 
the special case of quasistationary emission along an 
I ~ \m\ = 2 mode aligned with z, the eigenspaces of 
(^Li^g^Lb)) ^ are are z (eigenvalue to^ = 4) and the x — y 
plane (eigenvalue -I- 1) — rn^) /2 = 1). As quadrupo- 
lar emission dominates, these statements are a good ap- 
proximation for all time-domain averages (L(ciLb))c 

Small errors in high-order modes can be amplified by 
the factor particularly during the merger phase where 
higher harmonics become significant both physically and 
through numerical error. However, since higher-order I 
subspaces can be invariantly removed from ■04, a similar 
average can also be constructed for any subset of L To 
the extent that the dominant emission orientation is well- 
resolved, the directions derived from restricted averages 
should agree. 



III. OPTIMAL ORIENTATION VERSUS TIME 

The dominant principal axis of (L^qL^)^^ invari- 
antly determines a preferred orientation Z{t) versus 
time. Moreover, by design, any waveform that is well- 
approximated by quasistationary emission of a Z = |m| = 
2 mode along Z[t) must have Z{t) correspond to that 
instantaneous emission symmetry direction. On the con- 
trary, though a priori the eigensystem of the symmetric 
tensor (^L^^^Lb)) provides a frame, empirically the direc- 
tions associated with the two smaller eigenvalues can- 
not be trusted. During the inspiral the perpendicular 
directions are degenerate; small numerical asymmetries 
in each timestep cause the two smaller eigenvectors to 
jump chaotically in the plane perpendicular to the dom- 
inant eigenvector. Given the empirical degeneracy of 
the smaller two eigenvalues, we adopt a preferred frame 
at time i by a fiducial rotation: starting with a carte- 
sian frame aligned with the initial simulation's coordi- 
nate frame, rotate by i? = Rz{(l))Ry{0), which takes z 
to Z if 6*, are the spherical polar coordinates for Z. In 
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FIG. 1: Recovering a synthetic rotation: Starting with a 
equal-mass nonspinning (a — 0.0) binary waveform (bottom 
panel: "04;^ for I < 2, extracted at r = 60M), we apply 
a rotation taking z to the polar angles {6,(f>), where 6{t) — 
Ooil + 0.1cos27rf/P) and ^{t) = 2nt/P for P = 80M and 
9o = 7r/20. The bottom panel shows i/'4im after (dashed) 
and before (solid) this rotation, for the modes I = 2 and 
m = 2 (red), 1 (blue), and (green). Using (^L(^aLb))^, we 
recover an orientation Z(t) from the rotated waveform that 
agrees with the imposed rotation to within 5 x 10~® deg (top 
panel: 9{Z{t)) — 0{t)). Applying this rotation recovers the 
unperturbed waveform. A vertical dotted line indicates the 
peak I = |m| = 2 emission in both figures. 

terms of this rotation, the waveform in the instantaneous 
"aligned" frame at each time can be expressed in terms of 
the instantaneously "aligned" cigcnstates Y^~^\R~^n) 
and the simulation-frame expansion coefficients i^alm us- 
ing the representation theory of SU(2) [l3| : 

= Y.^ALMYtM\R-^n) (5) 

LM 

i^ALM = e^'^^''^^Y.^MApiR)^iLM' (6) 
M' 

The overall phase x does not enter into any calculation 
we perform and will be ignored. 

As a concrete example, in Figure [T] we compare the 
difference between a generic time-dependent rotation 
and the reconstructed orientation, obtained by apply- 
ing that rotation to a nonprecessing black hole binary 
simulation and reconstructing the optimal direction vcr- 
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FIG. 2: Aligning a processing binary: Starting with a 
precessing equal mass binary (ai — O.62, 02 = 0.6(3; 4- z)/V2; 
top panel), we derive the principal axis orientation from 
(!/((, ii,) restricted to / = 2 modes. The bottom panel shows 
the transformed harmonics. The orientation transformation 
is derived and applied at all times shown. To illustrate this 
method's robustness, in this figure we retain initial transients 
and late-time errors. 



sus time. The quasicircular orientation is recovered to 
(56* ~ 10~^ deg (i.e., 66^ comparable to working machine 
precision). As another example, in Figure [21 we compare 
the mode amplitudes for a spinning, precessing binary 
in an asymptotic incrtial frame (top) and in the "coro- 
tating" frame implied by this transformation (bottom). 
The transformed modes resemble the modes of a spin- 
aligned binary (Fig [TJ right panel). Note that in going 
from the generic-spin to "corotating" frame, the m ~ 1 
and m — Q modes have been reduced by more than an 
order of magnitude during the inspiral. 

Schmidt et al. Q have previously reconstructed opti- 
mal orientations from the emitted waveform. By contrast 
to their maximization-based method for tuning our 
algebraic method is fast, accurate, and invariant. This 
method can also be applied to any constant-^ subspace. 
For example, using the GT/PSU equal- mass, generic s pin 
simulation set summarized by O'Shaughnessy et al. |14l | 
we have reconstructed (i(ai(,)) for alH < 4 modes; for 
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I = 2 only; and for / = 4 only.^ All agree during the 
inspiral and at least some of the merger.'^ 

Finally, as this method operates on each I subspace as 
well as independently and algebraically at each timestep, 
it provides valuable, fast, invariant diagnostics of the 
quality of numerical waveforms. For example, numerical 
simulations can evaluate their waveform quality by the 
following three constructs: (a) eigenvalues of (^L(^aLb))^ 
should evolve smoothly and remain near their predicted 
values (i.e., near 4, 1, 1 if the I = 2 subspace is included); 
(b) the recovered dominant eigendirection of (^Li^^Li^-^^ ^ 
should be smooth in time; and (c) phcnomenologically 
speaking, the recovered paths from different multipole 
orders should agree. 



IV. OPTIMAL ORIENTATION VERSUS MASS 

In many cases of astrophysical interest, the emission 
direction changes significantly during the merger. As an 
example, BH-NS binaries can have their dominant emis- 
sion direction, the orbital angular momentum direction 
L, precess one to several times in a cone around the total 
angular momentum [lo| . In this case in particular and 
for high mass ratio binaries in general, the instantaneous 
emission direction averages out. Instead, the natural di- 
rection encoded by the signal is the total angular mo- 
mentum. Generally, however, a suitable average is not 
intuitively self-evident. 

A detector with comparable and coherent sensitivity to 
both gravitational wave polarizations is naturally char- 
acterized with a complex inner product that coherently 
accounts for both polarizations [Tsj . Specifically, two 
gravitational waveforms A, B expressed as curvature ('04) 
along a particular line of sight h are naturally compared 
by a complex inner product 



{A,B) 



df 



MfTBif) 



(7) 



Taking derivatives with respect to angle and expanding 
out the generators of the rotation matrix, the Fisher ma- 
trix {Tab = I {daip4:, dbipi) |/|(V'4, ■04)1) rcduccs to 



Tab 



|(£(aV^4,-^b)V'4) 
l(V'4,'04)| 



(8) 



The tensor l^L^^aLb)) therefore corresponds to the 



^ We have also tested the behavior of (L^^L^))^ using only the 
I = 3 mode subspace. During inspiral, this tensor has a well- 
defined dominant principal axis that agrees with the 1 = 2 re- 
sult, tlowever, at and beyond merger, (i(aii,))^ becomes nearly 
isotropic and does not provide a natural emission direction. 

^ A few tens of M after the merger, owing to their low and rapidly 
decaying amplitude, the I = 4 modes previously available are not 
resolved to the precision needed for directional reconstruction. 



orientation-averaged Fisher matrix.^ In other words, 
l^L^^aLb)) characterizes how well we can determine a 
binary's orientation, averaged over all possible directions 
along which we could see it. 

Moreover, for quasistationary precession dominated by 
a single mode at frequency /(<), the signal- weighted av- 
erage l^Li^aLb)) has a natural physical interpretation as 
a power- weighted average of (^Lf^aLb)) ^- Substituting the 
stationary phase approximation to the Fourier transform 



Mf) 



VWJdt 



(9) 



into the overlap {La4'ii Lbtpi) leads to an amplitude- 
weighted integral over (^L(^aLb))^ 



2dt 



-{LaMt)y{LbHm 



{2^f{t)YSH{f{t)r 

The orientation average of the Fisher matrix is therefore 



Tabdn/{^^) 



1 



l(V'4,04) 



2dt 



{L[aLb))t 



{2^f{t)YSHU{i 



■1) 



In other words, for sufficiently steady quasicircular pre- 
cession, <^i(aLb))jy/ is precisely the power-weighted aver- 
age of (L((jif,))j. The principal eigenvector of (L(aib))j^^ 
therefore corresponds to our intuition about an "average 
frame" . 



V. CONCLUSIONS 

In terms of orientation-averaged tensors, we have in- 
troduced an efficient algebraic method to extract fiducial 
orientations from waveforms extracted to infinity. Specif- 
ically, we have demonstrated that one natural orienta- 
tion of the emitted beampattern is associated with the 
dominant principal axis of l^Li^aLb))- The average can 
be calculated at any time or averaged, at any mass. In 
the time domain, our method accurately, efficiently, and 
smoothly reconstructs the emission direction throughout 
the inspiral. As with Q, if we retroactively align the 
waveform with this optimal direction at each timestep, 
the dominant modes become smoother and the subdom- 
inant modes significantly smaller. 

Conversely, if applied at each mass, our method also 
has close, provable connections to data analysis: it corre- 
sponds to the Fisher matrix for orientation angles, aver- 
aged over all possible lines of sight. In the case of extreme 
mass ratio, except for rare high-symmetry nonprecessing 
inclined orbits, our preferred direction reduces to the to- 
tal angular momentum, about which the orbit precesses. 



For brevity, we do not discuss maximization over coalescence 
time and phase; see [TBI for a more extensive discussion of the 
Fisher matrix derived from a complex overlap. 
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This invariant approach to the preferred orientation will 
benefit search strategies for numerical relativity wave- 
forms from strongly precessing binaries that coherently 
emplo y m ultiple modes, analagous to the proposal by Pan 
et al. for strongly precessing BH-NS binaries. 

We have proposed one of many possible definitions for 
the "instantaneous emission direction." Of course, other 
plausible generating tensors Q can be constructed with 
La alone, each corresponding to a different way of weight- 
ing the mode amplitudes i^Aim- ^^i^l more can be added 
by broadening the space of operators. That said, our def- 
inition offers a so-far unique feature: a provable connec- 
tion to an astrophysically relevant quantity (the Fisher 
matrix) . 

Other astrophysically significant features of numerical 
relativity waveforms could be employed to determine pre- 
ferred directions. For example, a "generalized equato- 
rial plane" could consist of the set of directions along 
which gravitational radiation is locally linearly polarized 
(e.g., as the surface 9tln'04/-04* = 0). Like the aver- 
age (^L(^al^b)}j these directions can be defined in the time, 
frequency, and mass domain. These emission directions 
lead to parameter estimation degeneracies: equal amount 
of left- and right-handed radiation imply two possible bi- 
nary orientations are consistent with that signal. Though 
intriguing and invariant, these symmetry directions de- 
pend sensitively on delicate cancellations of all available 
modes. Less invariant but equally significant is a direct 
comparison ( "overlap" ) of waveforms emitted in all pos- 
sible directions, as would be seen by a gravitational wave 
detector. This method will be discussed in a subsequent 
paper [isj . in the broader context of generic waveforms 



from spinning binaries. 

Finally and as discussed previously by Schmidt et al. 
3, both radiation and merger physics are more easily 
modeled in an "aligned" frame. For example, at present, 
hybrid waveforms have been constructed primarily for 
spin- aligned binaries for a handful of harmonics fig . [l7| 
(cf. Sturani et al. [iBl)- By tabulating and model- 
ing the "aligned" -frame waveforms and the corotating 
frame itself, hybrids can be constructed for generic pre- 
cessing waveforms. As another example, previous stud- 
ies of nonspinning and aligned-spin binaries suggest that 
their low-order modes evolve in phase with each other 
(arg'04;^ cx TOarg?/'422) through inspiral and merger 
[ig . l20j . By analogy, when expressed in a an "aligned" 
frame, low-order modes of generic precessing binaries also 
seem to evolve in phase. We will address these and other 
physical properties of "aligned" -frame waveforms in a 
subsequent publication. 
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